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Using the d-p model, we demonstrate that the pseudogap, which is induced by the supercon- 
ducting fluctuation, plays key roles in the determination of the phase diagram observed in high-T c 
superconducting materials. We take the pairing interaction mediated by the spin fluctuation and 
calculate the superconducting transition temperature T c , the NMR relaxation rate 1/Ti and the 
single-particle spectrum by treating both the superconducting fluctuation and spin fluctuation 
in a consistent fashion. As temperature decreases, 1/T\T increases at high temperatures, and 
it reaches a maximum followed by a sharp drop in the underdoped region, due to the evolution 
of the pseudogap in the single-particle spectrum. The evolution is also consistent with those of 
ARPES experiments. 

KEYWORDS: pseudogap, NMR relaxation rate, superconductivity, superconducting fluctuation, spin fluctuation, 
superexchange interaction, d-p model, 1/Af-expansion, slave-boson technique 
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In the high-T c superconductors, various anomalous 
properties have been observed at higher temperatures 
than the superconducting transition temperature T c , es- 
pecially in the underdoped region. Key issues on these 
anomalous properties must be on the existence of the 
pseudogap Wrhich. r^ft^WWh suggested based on various 
experimen^.ErBIMBBQ'Qu) In-plane resistivities pi 1 Hall 
coefficients^ and static uniform spin susceptibilities^ 
show peculiar temperature dependences below a certain 
characteristic temperature Tn,B> where the pseudogap is 
believed to become appreciable. Recently, photoemis- 
sion experiments have revealed that the pseudogap is of 
the same size and has. the same k-dependence as the su- 
perconducting gap.Q'tP More recently, the temperature 
evolution of the pseudogap has been observed by the 
angle-resolved photoemission spectroscopy (ARPES). □* 

There are several different views regarding the origin 
of the pseudogap. For example, it may be interpreted as 
a spin-excitation gap .in the singlet RVB states obtained 
in the t-J model,t£rO orJ±_may be due to the Bose con- 



girBB 



or due to superconducting 



densation of_bc 
fluctuation^' Bl^ 

In the present study, we demonstrate, using the d-p 
model, that the pseudogap induced by the superconduct- 
ing fluctuation (SC fluctuation) plays key roles in the 
determination of the phase diagram observed in high- 
Tc superconducting materials. The pairing interaction is 
derived from the antiferromagnetic spin fluctuation (AF 
fluctuation) induced by the superexchange interaction 
J s , while the AF fluctuation is suppressed by the pseu- 
dogap in the single-particle spectrum induced by the SC 
fluctuation. Thus, we must treat the effects of both the 
SC fluctuation and the AF fluctuation in a consistent 
fashion. Thus, we calculate the superconducting transi- 
tion temperature T c , the NMR relaxation rate 1/Ti and 
the single-particle spectrum. T c has a maximum at the 



hole-doping rate 8 = 0.11. At high temperatures, \jT{F 
increases as T decreases, due to the development of the 
AF fluctuation. We assign To as a temperature at which 
the effects of the AF fluctuation begin to be appreciable. 
Then this Tq corresponds to that [determined by experi- 
ments on the Hall coefficient Rft.w As the temperature 
decreases further, 1/T.T reaches a maximum at a cer- 
tain temperature, denoted by T sg , for 8 < 0.11, followed 
by a sharp drop at T < T sg . These features, are consis- 
tent with the results of NMR experiments.tH 2 ]) The plots 
of these three characteristic temperatures, T c , To and 
T sg , as functions of 8, account for essential features of 
the phase diagram of the high-T c superconducting ma- 
terials. The single-particle spectral weights at the Fermi 
level decrease near the [n, 0] and [0, it] directions as T 
decreases at T < To, and eventually the Fermi surface 
survives only near the [w, n] direction at T < T sg . These 
features.are consistent with the results of ARPES exper- 
iments, u 

First, we describe our model and formulation. We take 
the simplest version of the d-p model as 



H = e 



N 

m=1 k m 



, c k„ 



N 
m—l 



N 



n; 



m=1 k m ,i 



+ h.c.}, (1) 



which is treated within the physical subspace where local 
constraints 

N 

Qi EE 4m d im + bfh = 1 (» = 1, 2, • • • , N L ) (2) 

m—l 

hold. In the above, m , di m and hi are annihila- 
tion operators for a p-hole, a d-hole and a slave bo- 
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son, respectively, and tj(k) = tj i exp(— ik • Rj), where 
i^. = 2£[1 — i(cosfc x a + cos k y a)]i and a is the lattice 
constant. The suffix m represents the spin-orbital de- 
generacy (m = 1, 2, • • ■ , AT) introduced by dividing the 
phase space for k into N/2 subspaces for k m . keeping the 
total degrees of freedom for spins and orbitals unchanged 
as k = S m k • F° r numerical calculations in the 
present study, we set N = 2. The local constraints in eq. 
(2) are strictly held when rae .calculate an expectation 
value of a physical quantityli§Ej 

A set of self-consistent equations for single-particle 
Green's functions of the leading order in the 1/N- 
expansion_was solved to yield the p-hole Green's function 
given by c3 Ej 



G(k,ie„) 



E 

7=± 
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with 







4 s 





+ 0JQ+ 7v/( £ p - ^o) 2 + 46 t k 2 ^ 



(3) 



(4) 



(5) 



where ujq and & are the energy and the residue of the 
pole in the slave-boson Green's function. The chemical 
potential fi is determined by 

1 



1 



E 

ko-7 



(6) 



where n and 8 are the total hole number and the doped- 
hole number per unit cell, respectively, and f(E) = 
[ exp{E/kBT} + 1 ] — 1 . Here, it is noted that the so- 
lution E^ denotes the in-gap state emerging inside the 
charge transfer gap, A = e p — Sd, upon doping carriers to 
a Mott insulator. The band width of the in-gap states is 
proportional to 6 in the underdoped reigon. The effects 
of the interactions among the in-gap Sjtates ^re of highpr- 
order terms in the l/A^-expansionSa'SESED'Ea'Ej) 
Recently, it has been shown that quasi-particle interac- 
tions via the superexchange-interaction play dominant 
roles in underdoped systems & 

Now, we derive coupled equations to treat the effects 
of the SC fluctuation and the the AF fluctuation in a 
consistent fashion (termed the extended self-consistent 
i-matrix approximation). The pairing interaction medi- 
ated by the AF fluctuation via the superexchange inter- 
action J s in RPA is given byEJ 



V(q) = - 

1-J s (q)^ 0) (q,0) 

with J s (q) = -J s (cos(q x ) +cos(q y )), 



(7) 



Xi 0) (q,c + i0+) 
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It was shown that a component with the d x 2_ y 2 symme- 
try among various components of the spin-fluctuation- 
mediated interaction contributes dominantly to the pair- 
ing interaction^ Therefore, in the present study, we 
take only the component with the d x i_ y i symmetry of 
the pairing interaction as 



Vd 



-l^^^(k)^(k-k')^(k'), 
k k< 



(9) 



with ipd{k) — cos(k x a) — cos(k y a). 

The pairing susceptibility x(q> w )> which is the irre- 
ducible part of the i-matrix, is given by 

u2j-2 j-2 

1 k k 1 f 

x(q ' w + i0+) = -n- l ^ Mk)2 wn*f^ «L dx 

x [/(af)ImG(k, x + iO+)G(q - k, u - x + i0+) 

-/(-x)G(k, oj-x + iO+)ImG(q - k, x + iO^) 

A self-energy correction due to the SC fluctuation in the 
i-matrix approximation is given by 

-i Z-i Z i i />oo 

S(k, „ + i0+) - -^M*Y (( ^ - J__ dx 

x [f(x)X(q, x + uj + iO+)ImG(q - k, x + i0+) 
-g{x)lmX(q, x + iO+)G(q - k, x - u - iO^l)}) 

with 

X(q, uj + i0+) = ^ X (q, w + iO+)/(l - v dX (q, uj + 10+)). 

(12) 

The renormalized Green's function is given by 
G(k,cj+iO+) = [G{k,uj+iO + )- 1 -Y,{k,uj+iO + )}- 1 . (13) 



Now we replace G in eqs. (10) and (11) by G given 
in eq. (13). Then, eqs. (10)-(13) form a set of self- 
consistent equations for a given pairing interaction Vd- 
This procedure is called the self-consistent i-matrix ap- 
proximation. However, it is insufficient for the present 
case. The reasons are as follows. The single-particle 
spectrum obtained above has a pseudogap, i.e., suppres- 
sion of the density of states near the Fermi energy in the 
normal state. The pseudogap suppresses the AF fluc- 
tuation and hence the pairing interaction. Thus, these 
effects must be treated in a fully consistent manner. To 
this end, we take the following scenario. We calculate 
Xs (li 0) m ec l- (^) where G is replaced by G obtained 
above. Then, we find the pairing interaction Vd modified 
by the pseudogap through eqs. (7) and (9). We re- 
peat the entire process until we obtain a self-consistent 
solution for v d . This procedure is called the extended 
self-consistent i-matrix approximation. 

Using the solutions obtained above, we determine the 
_superconducting temperature T c by the Thoulcss crite- 



v k+q 



2 (E^ - w ) 2 rion, 



dxf(x) x [ImG(k + q, x 
+G(k + q,w + 



iO+)G(k,-^ + 
+ iO+)ImG(k, x + i0\ 



i0+) 



v dX (0,0)-v dX h(0,0) = 0, 



(14) 



where Xh is due to higher-order terms of the SC fluc- 
tuation, which is important in high-T c superconducting 
materials because of the quasi 2-dimensionality of Cu02 
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electronic systems. In the normal state, close to the su- 
perconducting phase, %h(q, 0) has a sharp q-dependence, 
practically, UdXh(qiO) is of 0(1) only at q = and is 
zero otherwise in the limit of T — > T c . Thus, \h does 
not contribute to the q-summation in the self-energy E 
(eq.(ll)) in this limit. In the present study, we simply 
assume that VdXh(0,0) is 0.5. 

We also calculate the NMR relaxation rate, 1/Ti, using 
the formulaH 

= F^^iSJ Im *q PA MM (15) 

with 

Xcf» = xi 0) (q^)/[l - Js(q)xi 0) (q^)], (16) 
-Fab(q) = [A a b + 2B(cos q x a + cosq y a)} 2 /2 

+ {A c + 2B(cosq x a + cosq y a)} 2 /2, (17) 

A ab = -170, A c = 20 and B = 40. (18) 

Here, Xs^(q, is given by eq. (8), where G is replaced 
by G. 

In actual numerical calculation, throughout the 
present study, we set 2t = 1.0 (which is of O(leV) in 
real systems), A = 2.5 and J s = 0.1. The total number 
of discrete points taken for the q-summation over the 2D 
first Brillouin zone is 32 x 32. The ui- integral over the 
region from — 2wo to 2u>q is replaced by the u-summation 
of 80 discrete points. 

Numerical results for the T-dependence in \jT\T are 
shown in Fig. 1 for the hole-doping rate S — 0.06 ~ 0.19. 
At high temperatures, 1/T]T increases as T or 6 de- 
creases, due to the development of the AF fluctuation. 
We define T as the characteristic temperature where the 
effects of the AF fluctuation begin to be appreciable. Ex- 
plicit values of To are determined by the condition that 
[&(l/T 1 T)\ T ^r ]/[{l/T 1 T)\ T ->oo,6=o.i] = -2 x 10 2 . At 
S < 0.11, 1/T]T has a maximum at a certain temper- 
ature denoted by T sg . Curves of \/T\T on the lower 
temperature sides terminate at T c for 8 > 0.1, but we 
have not yet obtained T c for 5 < 0.09 due to the poor 
convergency of the numerical calculation at T = T c . In 
the inset, the T-dependence of 1/TiT obtained at 5 = 0.1 
as above (solid line) is compared with that calculated by 
neglecting the effects of the SC fluctuation (dashed line). 
We note that the SC fluctuation also becomes apprecia- 
ble at T ~ To and continues to develop along with the 
AF fluctuation as T decreases, while the SC fluctuation 
induces the pseudogap in the single-particle spectra (see 
Figs. 3 and 4), which suppresses the AF fluctuation. 
Eventually, the developed pseudogap eliminates the AF 
fluctuation, which leads to a sharp drop in 1/TiT at 
T 5; T sg . Recently, Kontani showed that the Hall co- 
efficient i?H is also enhanced by the AF fluctuation. M> 
Thus, the initial enhancement of at high tempera- 
tures can be assumed to have the same origin as that of 
1/T\T. In this sense, we may argue that To obtained in 
the present study corresponds to T determined by i?H 
experiments.!]'' 

We plot the three characteristic temperatures To, T sg 
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Fig. 1. T-dependences of 1/TiT obtained with J a = 0.1 and A = 
2.5 for the hole-doping rate S = 0.06 ~ 0.19. In the inset, l/T\T 
(solid line) is compared with that calculated by neglecting the 
effects of the SC fluctuation (dashed line), where 8 = 0.1. 



and T c as functions of 6 in Fig. 2. Note that T sg ap- 
pears at S = 0.11 and increases as 5 decreases, whereas 
To appears at S = 0.26 and is much higher than T sg . 
We also note that T c has a maximum at S = 0.11. The 
dotted curve of T c at 5 < 0.09 is determined by eq. (14) 
with the help of the extrapolation of numerical results for 
the T-dependence of (1 — Vdx(0, 0) — VdXh{0, 0)), because 
we had difficulty in obtaining fully convergent solutions 
for the self-consistent equations (7)-(13) at T = T c for 
6 ^ 0.09. As S decreases, the quasi-particle band narrows 
and the nesting effect increases, which enhances the AF 
fluctuation and then pushes To higher. The AF fluctu- 
ation enhances the pairing interaction and then the SC 
fluctuation, which develops the pseudogap and then sup- 
presses T c in the underdoped region. T sg is determined 
by the competition between the AF fluctuation and the 
SC fluctuation. The phase diagram shown in Fig. 2 ac- 
counts for essential features of the anomalous metallic 
phase observed in high-T c superconducting materials & 

The single-particle density of states, p(u>) = 

-^mi^Ek^C 4 '^ + i0+ )' at T = 0.02, 0.025, 0.03 
and 0.2 with (5 = 0.1 where To = 0.057 are shown in Fig. 
3. At T = 0.2, p(u>) is nearly equal to that of the bare 
in-gap states without the effects of the SC fluctuation. 
At T < T , p(tjj) decreases near the Fermi energy (cj = 0) 
and increases at \u\ Z 0.05 as T decreases. The size of 
the pseudogap is roughly 0.1, which corresponds to the 
strength of the pairing interaction of order of J s . 

Lastly, we show contour plots of the k-dependence on 
the first Brillouin zone of the single-particle spectrum at 
the Fermi energy p(k, 0) = — ilmG(k, 0) with i5 = 0.1 at 
T = 0.2 in Fig. 4 (left) and atT = 0.02 in Fig. 4 (right). 
The single-particle spectrum p(k, 0) has a large value in 
bright regions, where the energy of a quasi particle with 
k is close or equal to the Fermi energy. Thus, the bright 
regions represent the shape of the Fermi surface. In the 
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J s =0.1 A=2.5 





Fig. 4. Contour plots of k-dependences (first Brillouin zone) of 
the single-particle spectra at Fermi energy p(k, u) = 0) with S = 
O.i, J s = O.i and A = 2.5, where T = 0.2 (left) and T = 0.02 
(right). Bright regions represent the shape of the Fermi surface. 



Fig. 2. S-T phase diagram with J s = 0.1 and A = 2.5. 
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Fig. 3. Single-particle density of states at T = 0.02, 0.025, 0.03 
and 0.2 with 6 = 0.1, J B = 0.1 and A = 2.5, where T = 0.057. 



dark regions, p(k, 0) is nearly zero, where the energy of 
a quasi particle with k is far from the Fermi energy. The 
Fermi surface exists in all directions at T = 0.2. As the 
temperature decreases, the single-particle spectra near 
the [tt, 0] and [0, 7r] directions decrease due to the SC 
fluctuation with d x 2_ y 2 symmetry, and then the Fermi 
surface remains only near the [ir, tt] direction at T = 0.02. 
These features are consistent with the results of ARPES 
experiments.& ) 
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